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Abstract 

The Riemann walk is the lattice version of the Levy flight. For the 
one-dimensional Riemann walk of Levy exponent < a < 2 we study 
the statistics of the support, i.e. set of visited sites, after t steps. We 
consider a wide class of support related observables M(t), including 
the number S(t) of visited sites and the number I(t) of sequences of 
adjacent visited sites. For t — > oo we obtain the asymptotic power 
laws for the averages, variances, and correlations of these observables. 
Logarithmic correction factors appear for a = | and a = 1. Bulk and 
surface observables have different power laws for 1 < a < 2. Fluctu- 
ations are shown to be universal for | < a < 2. This means that in 
the limit t — > oo the deviations from average AM(t) = M(t) — M(t) 
are fully described either by a single M independent stochastic process 
(when | < a < 1) or by two such processes, one for the bulk and one 
for the surface observables (when 1 < a < 2). 
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1 Introduction 



The Levy flight is a random walk in continuous space whose step size dis- 
tribution has a power law tail and is therefore sometimes called a "Levy 
distribution". The ubiquity of such distributions has been emphasized by 
many authors, and is a consequence of the power law tail being invariant 
under convolution. Many interesting instances of the occurrence of Levy 
distributions are given by Tsallis Q and Tsallis et al. [||. These range from 
applications in physics (superdiffusion, chaotic fluid flow) and engineering 
(leaking taps) through studies of the physiology of heart activity, all the way 
to descriptions Q of fluctuations of financial markets. 

A one-dimensional lattice version of the Levy flight may be constructed 
as follows. Let a random walk consist of independent steps, and let the 
probability p{£) for a displacement of £ lattice units in a single step be given 
by p(0) = and 

p{l) = A\l\~ 1 - a (£ = ±1,±2,...) (1.1) 

Here a > is the Levy exponent. Normalization of p implies that A^ 1 = 
2£(1 + a) where £ is the Riemann zeta function. This random walk was 
first studied by Gillis and Weiss || in 1970. It is called the Riemann walk 
by Hughes (Ref. |J, p. 154) and we will conform to that terminology. More 
generally we call of Riemann type any one-dimensional lattice walk whose 
p(£) is asymptotically proportional to |£|~ 1_a when \£\ — » oo. 

Riemann type walks were reviewed in detail by Hughes [||. Of partic- 
ular interest is the exponent regime < a < 2, where these walks have a 
mean square displacement per step, (£ 2 ), which is infinite. There then exists, 
at least for certain global walk features, a correspondence between simple 
random walk on a d-dimensional lattice and one-dimensional Riemann type 
walks of exponent a = 2/d. In some ways the fraction ^ acts as the walk's 
effective dimensionality. But whereas analytical results for noninteger di- 
mension d cannot be checked by computer simulations, the full continuum 
of a values is accessible to Monte Carlo studies. 

Much interest has centered around the following question. Let there 
be a t step Riemann walk. Then what are the statistical properties of its 
support S(t), i.e., of the set of sites that the walk has visited? There appears 
immediately an important difference between the exponent regimes < a < 
1 and 1 < a < 2. In the former regime the Riemann walk is transient |g, 0] 
and it is easy to show (see Sec. |2.5| ) that <S(oo) is a set of fractal dimension 
ds = a. In the latter case the Riemann walk is recurrent M, S(oo) 
coincides with the full one-dimensional lattice, and ds = 1. The existing 
literature deals with the different question of finding the properties of S(t) 
for asymptotically t; the results reflect, nevertheless, the same distinction 
between < a < 1 and 1 < a < 2. The borderline case a = 1 is more 
subtle. 
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Gillis and Weiss ||] study the number S(t) of distinct sites in the support. 
They find, among other results, that for t — > oo the average of this random 
variable behaves || as S(t) ~ t for < a < 1 and as S(t) ~ t}l a for 
1 < a < 2, where ~ indicates asymptotic proportionality. For a = 1 and 
a = 2 power laws with logarithmic correction factors appear ||. For a > 2 
the result S(t) ~ i 1 / 2 is identical to that for the simple random walk in 
d = l. 

A recent extension of this work is due to Berkolaiko et al. Q. Pursuing 
a question initially asked for the case of the simple random walk by Larralde 
et al. ||], these authors investigate the number S^it) of distinct sites visited 
by N independent t step Riemann type walks all starting on the same lattice 
site. Again power laws appear, both for t — > oo at fixed N and for N — > oo 
at fixed t. 



The present work extends the investigations of Gillis and Weiss into 
a different direction. We limit ourselves to the Riemann walk defined by 
Eq. ( |1 . 1|) , with a in the regime of greatest interest, that is, < a < 2. Our 
results may be summarized under three headings. 

1. Variance AS 2 (t). For any quantity X(t) we will denote its instan- 
taneous deviation from average by AX(t) = X(t) — X(t). Traditionally in 
this field the calculation of the average number S(t) of distinct sites visited 
has been followed by a calculation of the variance AS 2 (t) of that number. 
Thus, for the simple random walk S(t) was first calculated by Dvoretzky 
and Erdos @ in 1951, and AS 2 (t) by Jain and Pruitt in 1970. For 
the one-dimensional Riemann walk the present work supplements the 1970 
results due to Gillis and Weiss Q for S(t) by the corresponding ones for the 
variances AS 2 (t) in the regime < a < 2. 

2. Variables other than S{t). We use the powerful generating function 
method (GFM), which was introduced into the field of random walks by 
Montroll (l^] and Montroll and Weiss |flS|] . Overviews of this method are 
given by Weiss 0] and by Hughes [|| . The first calculation of a variance by 
the GFM, viz. that of S(t) for the simple random walk, is due to Torney 
[b§ in 1986. 

In 1994 Coutinho et al. [15] performed Monte Carlo simulations of, 
among other things, the number of unvisited islands enclosed by the support 
of the t step simple random walk in two dimensions. This led Caser and 
Hilhorst [16] to analytically determine the asymptotic behavior of the aver- 
age number of islands. Subsequently Van Wijland et al. [17, 18] developed 
a compact GFM based analytical scheme for calculating simultaneously the 
averages, variances, and correlations of a large class of observables charac- 
teristic of the support, generically denoted by the symbol M(t). In d = 2 
this class includes also the total boundary length of the support, and in 
d = 3 its surface area and Euler index. 

Here we bring this scheme to bear on the one-dimensional Riemann walk. 
The support of this walk consists of alternating sequences of visited and un- 
visited sites. Among the most prominent members of the class of observables 
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M(t) is, next to S(t), the number of visited sequences, that we will denote 
by I(t). Islands in d = 1 just are unvisited sequences enclosed by the sup- 
port, of which there are I(t)—1; the support furthermore has 2I(t) boundary 
sites (= visited sites adjacent to an unvisited one). Table I summarizes our 
results for the asymptotic laws of the averages, variances, and correlations 
involving S(t) and I(t). Beyond their intrinsic interest these laws may serve 
in heuristic arguments in reaction-diffusion processes, e.g., to estimate the 
trapping probability of an atom that diffuses in a random absorbing envi- 
ronment, or the effective reaction rate between two diffusing species. We 
defer further comments to Sec. 0. 
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Table I. Leading asymptotic behavior as t — > oo of the averages, variances, and 
correlation of S(t) and I(t) in different regimes of the Levy exponent a. The exact 
prefactors of the asymptotic laws are given in the text. The results for S(t) are due 
to Gillis and Weiss ||; the result for AS 2 (t) in the range < a < | follows from 
the theorem of Jain, Orey, and Pruitt (see Hughes g, p. 344); all others are new. 



3. Universality of fluctuations. The deviations from average AS(t),AI(t), 
. . . ,AM(t),. . . are randomly time-dependent variables that one would a pri- 
ori expect to exhibit some degree of correlation. One calls these fluctuations 
universal - by lack of a better name - when in the limit t — > oo all AM(t) 
are asymptotically equal (up to a proportionality constant) to a single M 
independent stochastic process. For the simple random walk universality 



was shown to hold in dimensions d = 2 [17] and d = 3 [18], and not to hold 



in d = 4, 5, . . . . For the d = 1 Riemann walk we find that universality holds 
in the exponent regime | < a < 2, but not for < a < |. A novelty with 
respect to the case of the simple random walk is that for 1 < a < 2 not a 
single, but two M independent processes are needed to describe the univer- 
sal fluctuations: one applies to bulk and the other to surface observables. 
The precise statements are given in Sec. ||. 

This article is set up as follows. Sec. [2] describes those elements of our 
analysis that are common to the full exponent interval < a < 2. Secs.^ 
and H deal more in particular with the exponent regimes < a < 1 and 
1 < a < 2, respectively, and derive the asymptotic behavior of averages, 
variances, and correlations. In Sec. [5] we do the same for the exceptional 
values a = | and a = 1. In Sec.^ we discuss the universality properties. In 
Sec. [7| we provide some additional interpretation of our results and conclude. 
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2 Observables, averages, and correlations 



2.1 Observables M(t) 

Our analysis is based on first writing quantities of interest in terms of the 
field m(x, t) of "complementary occupation numbers" denned by m(x, t) = 1 
if site x has not yet been visited at time t, and m(x, t) = otherwise. The 
expressions of S and I in terms of m are 

oo 

S(t) = J2 (2-1) 

x=— oo 
oo 

I(t) = m(x,t)[l-m(x+l,t)] (2.2) 

30=— OO 

S and I are representatives of a general class of "observables" M that are 
sums on x of a summand to which each lattice site contributes a factor m 
1 — m, or 1, i.e., the summand tests for the presence of a specific pattern of 
visited ("black") and unvisited ("white") sites. The following slightly more 
abstract characterization of the M will be needed. Let A = {a} be a finite 
set of distinct nonnegative integers a, such that either A = or, if not, 
A includes the element a = 0. The general observable M(t) that we will 
consider is 

oo 

m w= e E^n^+M) ( 2 - 3 ) 

x=— oo A a€A 

where for A = the product is equal to unity and where {[ia] is a set of 
numerical coefficients characteristic of M. When their M dependence needs 



to be indicated we will write ha[M]. Eqs. fl2.2|) and fl2,l|) show that S(t) 
and I(t) are of the form of Eq. ( p.3|) with only two nonzero coefficients, as 
shown in Table II. 



A 


Ha[S] Ha[I] Ha[Si] Ha[Ii] 




{0} 
{0,1} 
{0,2} 
{0,1,2} 


1 

-1 1 1 
-1 -2 
1 

-1 1 



Table II. Coefficients /^[M] for the four observables M = S,I,Si,Ii denned in 
the text; entries not shown are zero. The coefficients in each column add up to 
zero. 

Two further examples of observables of type fl2,3|) are the total number 
S±(t) of visited sequences consisting of only a single site, and the total num- 
ber I\(t) of single-site unvisited sequences. Their coefficients ha involve sets 
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A of up to three elements; they are easily determined and have also been 
listed in Table II. 

The following remarks, important for later, are verified without much 
effort. The coefficient is nonzero if and only if M is built up exclusively 
out of factors m. Since these correspond to visited sites, that make up the 
"bulk" of the support, we will call an M of this type a bulk observable. 
Observables built up exclusively out of factors m do not occur, since their 
expectation value on an infinite lattice is infinite. Hence the remaining 
observables refer to patterns consisting of both visited and unvisited sites, 
and we will therefore call them surface observables. [In the terminology 



of Refs. [17, [j^] these are "black" and "black-and-white" observables. They 
might also be called "S-like" and "I-like", respectively] The distinction 
between these two subclasses will play a role only in the exponent regime 
1 < a < 2. 

2.2 Basic formulas for averages and correlations 



In this work we will first evaluate the t — > oo behavior of the averages M(t). 
Then we turn to the covariance matrix AM(t)AM'(t), where M'{t) is a 
second observable with coefficients fi' A . Although the authors of Refs. jT?]] 
and (l8j deal with the simple random walk, the larger part of their formal 
developments also holds for the Riemann walk. 



The averages M(t) and M(t)M'(t) can be obtained as follows [17, 



Let G(x, t) be the Green function of the one-dimensional Riemann walk, 
that is, the probability for a walker starting at the origin to occupy site x 
after t steps. Let G(x, z) = Ylu=o z t G(x, t) denote its generating function 
and let Ga(z) be the \ A\ x \ A\ matrix of elements G(a — a', z) with a, a' £ A. 
From this matrix one constructs the "inverse sum" Ga(z) defined by 

G A \z) = \G- A \z)U> (2-4) 

a,a'<=A 

These scalars satisfy certain elementary relations stated in Appendix A as 
Properties 1-3. Two functions Cm(z) and Cmm'(z) are defined in terms 
of the Ga(z) according to 

C M (z) = £ (2.5) 

00 1 11 
Cmm'(z) = YjYjUAh'b J2 \r FT ~ r f \ ~ r r \ ] ( 2 " 6 ) 



Here A U (r + B) denotes the union of the set B, translated by r, and A. 
The averages M(t) and M(t)M'(t) are then obtained as @ 



where the integrations are counterclockwise around the origin. 



The identities (2.4)— (2.8) are fundamental to random walk theory; they 
hold for any translationally invariant random walk, whether on a finite lat- 
tice with periodic boundary conditions or on an infinite lattice. They allow 
for the calculation, in a very compact way, of many known and new results. 

Special cases. When M(t) = S(t), the following simplifications occur. 
The sums on A and on B in Eqs. ( |2.5| ) and ( |2.6| ) then have only the single 
term with A = {0} and B = {0}, respectively. Furthermore G{o}(z) = 
G(0,z), the matrix G{ }u(r+{o}) ( z ) 1S two by two, and an easy calculation 
leads to G{ }u(r+{o}) (^O = |(G(0, z) +G(r,z)). When M(t) = I(t), the sum 
in Eq. ( fT5| ) involves GjmQ?) and G{ 0) i}(z) = ^(G(0, z) + G(l, z)). The sums 
on A and B in Eq. ( |2.6| ) then lead to four terms, which may be evaluated 
with a little more effort. 



2.3 Limit t — > oo and scaling limit 

Explicit evaluation of the general expressions Q2,4| ) -( |2.8| ) is limited in prac- 
tice by the calculation of the inverse sums G^, which require the inversion of 
a matrix of dimension \A\. Similarly, evaluation of Ga\j(t+B) is an inversion 
problem of dimension \A\ + \B\ (when A and r + B have an empty intersec- 
tion). It turns out that the sum on r in Eq. (|2.6| ) can be performed only in 
the scaling limit 

z — ► 1, |r| — > oo with £ = r(l — z)« fixed (2-9) 

Finally, it will be possible to evaluate the integrals in Eqs. ( p. 7]) a nd ( |2.8| ) 
only asymptotically for £ — > oo, a limit already implied by Eq. (|2.9|) . 

In order to prepare for these limits we rewrite the preceding expres- 
sions as follows. Using the simplified notation Gq(z) = G(0,z) we split the 
generating function G(x, z) up according to 

G(x,z)=G (z)-g(x,z) (2.10) 

In full analogy to Ga(z) we define Sa( z ) as the matrix of elements g(a—a', z) 
with a, a' G A, and g^ 1 as the sum of all elements of g^ 1 - Let now J be the 
square matrix of elements J aa i = 1. Then 

G A (z) = Go(z)3 - Sa (z) (2.11) 

and, by Property 1 of Appendix A, 

G A (z) = G (z) - g A (z) (2.12) 

This splitup will be useful for studying the z — > 1 behavior of Ga(z). Al- 
though Gq(z) may (1 < a < 2) or may not (0 < a < 1) diverge as z — > 1, 
the functions g(x,z) and gA( z ) remain finite in that limit. 
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We now turn to the inverse sum G A{J ^ t+ b) constructed from the matrix 
Gj4u(r+B)' The dimension of this matrix is typically \A\ + \B\. Let 3 AB be 
the \A\ x |5| matrix with all J^ b B = 1. We then have (for i fl fi = 0) 



/ G A (z) 6(r,z)J^ \ / V 
GAUir+B){z) -{G(r,z)3^ G B {z) ) + {v T 



(2.13) 



where V is the matrix of elements 

V a , r +b = G(r + b-a,z)- G{r, z) a £ A, b e B (2.14) 
and \~ T is its transpose. The first matrix on the RHS of ( 2.1 3| ) has the form 



( |A.4[ ) of Appendix A. Applying Property 3 to that matrix we conclude 
that 

1 _ G A (z) + G B (z)-2G(r,z) + 2 



&Avj{r+B){z) G A {z)G B (z) - G 2 (r,z) ' / 2 _ 15 ) 

where we anticipate, and will have to show later, that V is small, that the 
correction terms are of order C(V 2 ), and that they are negligible for our 
purpose. 

Further analysis depends on the exponent a. We consider the two main 
regimes < a < 1 and 1 < a < 2 in Sees. [3] and ||, respectively. The excep- 
tional values a = | and a = 1 are discussed in Sec. ||. 



2.4 Riemann walk Green function 

All quantities of interest have been expressed above in terms of the Riemann 
walk Green function G(x,z). An elementary calculation yields 

f w da e ~ iqx 

G(x,z) = / — w r (2.16) 

-i oo 

The q — s- behavior of X(q) is crucial for the large scale features of the 
Riemann walk. It is known j|, |(| that 

X(q) = 1-C a \q\ a + 0(q 2 ) (q -» 0) (2.18) 

where for completeness we state the explicit expression 

C~ l = 2C(1 + a) T(l +a)/[ir sin(avr/2)] (2.19) 

One finds by standard methods (see e.g. ||) that in the limit z — > 1 the 
Green function in the origin Gq(z) has the asymptotic expansion 

G (z) = G (l) - B a (l - z)^- 1 + 0(1 - z) (0<a<l;a^ \) (2.20) 

G (z) = ilogKl-zJ-^ + Oa-^) (a = l) (2.21) 
G (z) = A a (l -z)- 1+ i + O(l) (Ka<2) (2.22) 
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where B a and A a are the constants 

B a = -Cl /a /[2sm(ir/a)] (± < a < 1) (2.23) 
A a = \/[2aC l J a sm{^/a)] (1< a < 2) (2.24) 

and c is a constant such that there is no 0(1) term in Eq. ( p. 21 ). In 
Eqs. (|2T2^ )-(l2~2l) and elsewhere we use the following convention. The sym- 
bol 0(X) indicates terms that are of order X in the applicable limit 
(X — > or X — > oo); this however is not to say that all preceding terms 
are larger. Thus, the nonanalytic term in Eq. ( |2.20| ) is larger than the 
0(1 — z) terms only for ^ < a < 1. For < a < \ it is present only as 
a correction to the 0(1 — z) terms; the expression for its coefficient B a 
in that regime is different from Eq. ( |2.23 ) but will not be needed. For 
the borderline case a = ^, excluded from Eq. ( 2.20| ), we have Gq(z) = 
Go(l) — \ (1— z) log(l — z)^ 1 + 0(1 — z); but this special nonanalytic behavior 
will stay subdominant everywhere in the remainder. 

From Eqs. tfTWj ) and pl8| ) one deduces that in the scaling limit ( f2~9| ) 

G(r,z) ~ (l-z)"" 1 ^^) (0<a<2) (2.25) 

where ^ = r(l — z)& and F(^) is the scaling function 

For ^ — > it behaves as 

F(f) ~ 2a({l + a)/[TTsm(aTr)} C 1+a (0 < a < 1) (2.27) 
F(0 ~ ilogT 1 (a = l) (2.28) 
F(0 = A^ + O^ 1 ) (l<a<2) (2.29) 

In Sees. 4 and 5 we will also use the function f(r, z) defined by 

G(r, z) = G (z)f(r, z) (1 < a < 2) (2.30) 

In the scaling limit one has f(r, z) ~ /(^) = F(£)/F(0) when 1 < a < 2. 



2.5 Support at £ = oo 

Whereas the remainder of this paper deals with the large t behavior, we 
briefly comment here on the structure of the support S(t) at t = oo. 

As is well-known J?], ^]), random walks are recurrent (are transient) if 
Gq(1) = oo (if Go(l) < oo). The Riemann walk of this work is recurrent 
for 1 < a < 2, which means that all sites are visited with probability 1, and 
that at t = oo the support «S(oo) coincides with the full one-dimensional 
lattice. 
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For < a < 1, however, the Riemann walk is transient, so that at 
t = oo the support «S(oo) will still be only a subset of the full lattice. We 
may estimate the average number of visited sites Sj, between x = —L and 
x = L in 5(oo). According to standard random walk theory [[?], |j 



Upon substituting ( p. 16 ) in ( 2.31| ) one easily evaluates £l for asymptotically 



large L, with the result that £_l ~ L a . It follows that the support S(oo) 
has fractal dimension ds = a. 

3 Riemann walk of exponent < a < 1 

3.1 Averages 



The large time behavior of M(t) comes from the behavior of Cj^(z), defined 
in Eq. (|2.5| ), in the limit z — ► 1. From Eq. ( [2.10D and the explicit expressions 
( |2.16| ) and Q2.17 ) it may be shown that g(x, z) = g(x, 1) + 0(1 — z) for all 



< a < 1. Hence 

<7a(z)=0a(1) + 0(1-*) (3-1) 
after which it follows from Eqs. ([DJ), (]2~2^ ), and (^TTJ) that 

G A (z) = G A (1) - B a (l - z)^~ l + 0(1 - z) (3.2) 
Inverting this relation and substituting in Eq. ( |2.5| ) gives 

C M (z) = -mi - B a m 2 (l - z)^~ l - B 2 a m 3 (l - z)^~ 2 + . . . + 0(1 - z) 

(3.3) 

where the m n are determined by the coefficients ha of the observable M 
according to 

m n [M] = E^y (n = l,2,... ; < a < 1) (3.4) 



In Eq. ([T^) the dots stand for a power series in (1 — z)q _1 and the number 
of nonanalytic terms between the zeroth and the first power of 1 — z is equal 
to n a = pj^jl — 1. That is, n a is zero for < a < \ and, as a goes up, 
jumps to 1, 2, 3, . . . at a = |, |, |, . . . , respectively. After Laplace inversion 
we get for the average M(t) in the limit t — > oo the explicit result 



M(t) = mi t + " m 2 t l -~ + % m 3 t^ + . . . + 0(1) 

V a/ I a/ (3.5) 

where the number of nonanalytic terms between the leading and the 0(1) 
term is again equal to n a . 
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3.2 Correlations 



In this subsection we consider two - possibly equal - observables M and 
M' , represented by sets of coefficients {ha} and {//4J, respectively, and 



wish to study their correlation. The starting point is Eq. (2.6) for Cmm'(z), 



in which we substitute Eq. (2.15). Whereas Ga(z) and Gb(z) tend to finite 



values in the limit z — ► 1, the Green function G(r,z) vanishes in that limit 
when taken with £ fixed. This suggests that we expand in powers of G(r, z), 



Gau{t+b){z) G a (z) G b (z) n=1 
with coefficients 



J2C^(z)G n (r,z) + 0(V 2 



(3.6) 



[G A (z) + Gb(z)][Ga(z)Gb(z)} ? 1 (n even) 

n + l 



-2[Ga(z)G b (z)}-^ (nodd) (3.7) 



The sum over space that occurs in Eq. ( |2.6| ) leads us to now consider the 
sums J2 r G n (r, z). From conservation of probability one finds that for n = 1 

Y^G(r,z) = (l-z)- 1 (3.8) 



For general n the calculation of J2r G n (r, z) is slightly more laborious; after 
substituting Eq. ( |2.16| ) for G one finds by explicit expansion in powers of 
1 — z that for z — ► 1 the sum on r behaves as 



G"(r, *) ~ F a , n (l - + O(l) (0<a<l; a ^ 1 - i) 



(3.9) 



For n = 1 — j^- (with n = 1,2,...) the sum on r instead behaves as 
log(l — z); this happens, in particular, for n = 3 when studied 
separately in Sec.||. 

We will now continue to consider the generic case. The nonanalytic term 
on the RHS of Eq. (|3.9|) dominates the 0(1) term only for n < j^— , In that 

case ( p.9| ) follows just from the scaling form ( |2.22 ) of G and from the £ — > 



behavior ( |2.27| ) of F(£). One then finds for the prefactor F a>n the expression 

POO 

F a , n = 2 / d£ F n (£) (1 - i < a < 1) (3.10) 



For < a < 1 — — the expression for F a ^ n is different and will not be 
needed. Eq. ( [3. 10 ) shows that the main contribution to the sum on r comes 



from £ ~ 1, that is, from r ~ (1 — z)~*. For n = 1 and n = 2 the integral 
(|3.10|) yields the explicit results F a A = 1, in agreement with Eq. ( |3.8D , and 
F a ,2 = (a ~ respectively. For n > •j^— the sum on r draws its main 

contribution from the short distance (nonscaling) regime r ~ 1, and is of 
0(1) for z -»■ 1. 
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By successively substituting Eq. ( j3?t ) in Eq. (3.6), neglecting the C(V 2 ) 
terms in that equation - which is justified in Appendix B - , then substitut- 
ing Eq. ( |3.6| ) in Eq. fl2.6|), expanding Ga{z) and Gb(z) according to Eq. (]37 
and using the z — > 1 behavior of J2 r G n (r, z) obtained in Eqs. (|3.8| ) and 
we find 

Cmm'(z) = -2(1 - z)~ 1 m 1 m 1 

- (1 - z)~~ 2 B a [3 - i](mim2 +771,27711) 



-(l-z)« 3 ^4 - ^JBl(mim' 3 + m 3 mi + m 2 m 2 ) + 2F a ^m 2 m' 2 
-... + 0(1) (3.11) 

Here the dots stand for terms of order (1 — z)~ l+k<y ~~ l \ with k = 3,4, . . . ; 
and the are related to M' in the same way as the m n are to M. For 
t —s- 00 we therefore find by substituting (3.11) in Eq. (|2.8|) and Laplace 
inverting 



1+2 



M(t)M'(t) = m\m' x t 



+ 



B 2 2F 
r(4- 2) ( mim 3 + m sK + m 2 m 2 ) + ^£ g 



a, 3 / 



4-2 

t a 



+ --- + 0(t) (3.12) 
The successive terms in the above series all have one power of t more than 



the corresponding terms in the series ( |3.5| ) for M(t), and the number of 
nonanalytic terms between the leading and the 0{t) term is once more equal 
to n a — 1. The product M(t) M'(t), which follows from Eq. (|3.5|), now has 
to be subtracted from the series ( |3.12| ). This exactly cancels the terms in 
(|3.12|) proportional to t 2 and to i 3 ~a but leaves those proportional to t 4 ~& 
and of 0{t). The f 4- ^ terms are leading only if | < a < 1. Hence we find 
for the correlation between observables M and M' in the limit t — > 00 



AM(t)AM'(t) ~ m2 m 2 t 4 ~ (3.13) 



valid for | < a < 1, and in which 



a r ( 5-^) + r(4-|) r»(3-i) ^ i4j 

We have supposed here that m 2 , m' 2 7^ 0. The preceding analysis changes 
when either of these two coefficients vanishes. We do not know of any 
physically interesting examples where this happens, and do not pursue our 
analysis in this direction. 
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The borderline case a = | is considered in Sec.||. In the interval < 
a < | the calculation of the present section applies, but with the result that 



AM (i)AM'(i) ~ K MM 't (0 < a < 



(3.15) 



in which the coefficient kmm' has contributions from the 0(1) terms in 
Eq. ( |3.5|) and the 0(t) terms in Eq. ( 3.12 ), and does not factor into an M 
and an M' dependent constant. This difference between Eqs. ( 3.13j) and 
( |3.15| ) is crucial for the phenomenon of universality discussed in Sec. ||. 



4 Riemann walk of exponent 1 < a < 2 

4.1 Averages 



The calculation of M(t) starts again from the series ( |2.5| ) for Cm( z )- The 
calculation in the exponent regime 1 < a < 2 is different from that of the 
preceding section because now Gq{z) diverges as z — > 1. Since gA (z) remains 
finite for z — > 1, this suggests that we use Eq. ( 2.12j ) and expand Ga(z) in 
powers oi gA(z)/Go(z). This yields 



Cm{z) = ^2 ha 



1 



9a(z) . 9 A ( Z ) g\ 



A ^ G (z) L 1 + G (z) + Gl{z) + ° { Gl - (41) 

We now substitute in Eq. Q4.1| ) the expansion ( 2.22j ) for Gq(z) and use that 
5a(1) is finite. The result is a power series in 1 — z in which there appear 
coefficients that we denote again by m n but that are defined for 1 < a < 2 
as 

m n [M] = -J2 »A9a(1) (n = 0,1,2,... ; 1 < a < 2) 

(4.2) 

It will turn out that we need only the leading term, which is 

^- 1 m (l - zf-T, (m + 0) 
C M {z) ~ <! 2 (4.3) 

A- 2 mi(l - zf~« (mo = 0,mi^0) 



We pause to note that in the terminology of Sec. 2.1 the condition mo 7^ 
characterizes the bulk or "S-Uke" observables, and the condition mo = 
the surface or "I-like" observables. This is the first equation where a 
difference appears between these two subclasses; in its analog, Eq. Q3.5| ) of 
the preceding section, no such distinction appears. 

Upon using Eq. (^4.3|) in Eq. (|2.7|) we obtain after an inverse Laplace trans- 



formation the asymptotic expansion of M(t) as t — > 00, 

f [A.ra + ^rw- (mo^o) 

M{t) ~ <^ (4.4) 
1 [AJX^-W" -1 (m = 0, mi ^0) 

where the dots indicate terms of lower order in t. 
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4.2 Correlations 



For the calculation of the correlation M(t)M'(t) via Eqs. (|2.8| ) and (|2.6| ) we 
have to return again to expression (|2.15| ) for ^/Q>a\j(t+b){z), which is needed 
in Eq. (|2.6|) . We use Eqs. ( 2.15|) , ( |2.12|) , and ( 2.3C| ) to rewrite this quantity 



MS 



9/1 f( r 7 \\ 9a(z) 9b(z) 

W-nr>z))-lgfa G^) +0{ y^ 



&Ali(r+B\{z) G (z) 1 _ f2( r z ) _ 9A(z) _ g B (z) g A (z)g B {z) > 

The function f(r,z) was defined in Eq. fl2.30p . In the scaling limit f(r,z), 
gA(z), and gB( z ) have finite limits, whereas Go(^) diverges. An expansion in 
inverse powers of Gq(z) corresponds therefore to an expansion in ascending 
powers of 1 — z. Writing for short /, gA, gs-, and Go when f(r,z), gA(z), 
gB(z), and Gq(z) are meant, we find after a straightforward calculation 



If r 

Cmm'(z) = E 2 +(2 + /) 

r=-oo 1 "T / A,B^0 



gA + gB 

Go 



2-2f-f g\+g 2 B 



Let us write = m n [M'\ for the coefficients that characterize the ob- 
servable M' . The two distinct cases described by Eq. ( |4.4j ) now lead to the 
following possibilities. 

Case (i): uiq 7^ and m 7^ 0. In this case the leading term in the 



expression in brackets in Eq. (4.6) survives under the sum on A and B. 

Case (ii): rriQ = 0, mi 7^ 0, and m' 7^ 0. In this case in order to survive 
a term in the bracketed expression should contain at least one factor gA (z) . 

Case (iii): m,Q = m' = but mi 7^ and m'i 7^ 0. In this case a term in 
order to survive must contain at least one factor gA(z) and one factor <7b(z)- 



Upon using in each of these cases for Go(z)the expansion ( 2.22j ), passing to 



the scaling limit, and writing / for /(£), we find that the result is 

-A~ l (l - 2) 1- «/oo m m' 
C M M>(z)~{ _A-i(l-2) 2 -|/ 10 mim (4.7) 
-A~ l (l - z) 3 -«/ u mimi 
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in the three cases (i), (ii), and (hi), respectively; here the coefficients fke 
represent the integrals 



foo 
ho 



f 
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/•CO 

4/ duii + fr 1 

Jo 

/•CO 

2/ dC/(2 + /)(l + /)- 1 
Jo 

/•CO 

4/ du(i-frHi + fY 

Jo 



(4.8) 



After substituting Eqs. (|4.7j ) in Eq. (|2.8| ) and carrying out the inverse Laplace 
transformation we find, in the limit t — > oo, 



M{t)M'(t) ~ { 



I T- 1 (-l + ^) A^f u m 1 m' 1 d- 2 



(4.9) 



respectively, for the three cases distinguished above. Upon combining these 
results with those of Section 4.1 one obtains, for t — > oo, 



' Bl 00 m m' t 2 c 



AM(t)AM'(t) - 

, B 2 all raxm! x t«~' 2 
in which the coefficients B a ke are given by 

Bloo = A'Hfoor-Hl + I) - A?T-*(1 + ±)] 



(4-10) 



J al0 

B 2 all = A-*[f n r-\-l + ±)- A^T~ 2 (j-)] 



(4-11) 



l-p-2/2' 



in the three cases (i), (ii), and (hi) defined above, respectively. We recall 
that the fke on the RHS of Eq. ( 4. 1 1| ) are given by Eq. ( |4.8| ) as integrals on 
/(£), with /(£) in turn given by Eq. (|2.26|) . 



Riemann walk of exponents a = \ and a 



5.1 Exponent a — 



3 



In this special case M(t) is still given by Eq. (3.5). However, the calculation 



of M(t)M'(t) has to be reconsidered, as signalled by the fact that F a ^ in 

2" 
3 



Eq. (3.11) diverges for a — > I .In order to calculate J2r G(r,z) we cannot 
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now use Eq. fl3.9|) . Instead we replace G(r, z) by its scaling form ( |2.22 ) but 



take into account that |£| has a lower cutoff |£| ~ est x (1 — z)-i. This gives 
YG 3 (r,z) ~ 2 / d£F 3 (0 

~ iC 2 log(l-z)^ 1 +0(1) (z->l) (5.1) 



where in the second step we used Eq. ( p. 24 ) and found for the coefficient 



the value C 2 = 2 12 3~ n / 2 7r~ 3 C 3 (§). In this case, due to the log(l - z) in 
the equation above, M(t)M'(t) is larger than the product M(t) M'(t) by a 
factor logi, and determines by itself alone the final result, which reads 



AM(t) AM' (t) ~ C 2 m 2 m' 2 t log t (t -»■ oo) (5.2) 
This tlogt behavior is the same as in the well-known case of the simple 



random walk in spatial dimension d = 3 [18]. 



5.2 Exponent a = 1 

The case of Levy exponent a = 1 is subtler than the others. Since it is 
closely analogous to the simple random walk in dimension d = 2 |17[] , we 
will not present all steps in detail. Eq. ( 2.21| ) shows that for a = 1 the 



Green function in the origin, Gq(z), diverges as z — > 1. We can therefore 



expand Cm{z) as a series in the same way as in Eq. (4.1). Since here again 
gA{z) = 5a(1) + 0(1 — z), and in view of the logarithmic behavior ( |2.21|) , 
this series now leads to an expansion of Cm(z) in inverse powers of Go(z). 
If the first nonzero term is of order k + 1, then we have explicitly 

C M {z) = -m k GQ k -\z) - m k+1 G k - 2 (z) - m k+2 G k - 3 (z) + ... 

(5-3) 

with the m n defined by Eq. (|4.2| ). The cases of physical interest have k = 
(bulk observables) or k = 1 (surface observables) , but it will be notationally 
convenient to keep A; as a parameter. We will also refer to it as the order of 
M. 

To find Cmm'( z ) we may still start from Eq. ( (4.5[) , but now the expansion 
of this equation runs differently. The reason is that for z — > 1 at fixed £ the 
function f(r,z) (defined by ( |2.30D ) behaves as F{£)/Gq(z) and so is of the 
same order as gA(z)/Go(z) and gB{z)/Go(z). We therefore have to perform 
a double expansion of the RHS of Eq. Q4.5P in terms of on the one hand 
F/Go and on the other hand qa/Gq and qb/Gq. The sum on r, which in 
the scaling limit becomes an integral on £, then leads to the appearance of 



coefficients F\ n defined as in Eq. ( 3.10 ) but with a = 1. Special cases are 
F\ t i = 1 and F\p = |. Let m k and m' k , be the first nonzero coefficients in the 
expansions of Cm{%) arid Cm'(z), respectively. Then we find for Cmm'( z ), 
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retaining only the three leading order terms in the limit z — > 1, 
CmM'( z ) — TTTTJTn 2m(,mL 

v ' (l- z )G k Q +k+2 { z y k 

-G 1 (z)(^(k + k' + 2)a 2 m k m' k , - 2(m k m' k , +1 + m k+ im' k ,)^ 

+G 2 (z)(2F 1 ^(k+l)(k' + l)a 3 m k m' k> 

— g(A; + k' + 3)a2(m k m' k , +1 + m k+ im' k ,) 

+ 2{m k m' k , +2 + m k+ im' k , +1 + m k+2 m k ,)) (5.4) 



The inverse Laplace transforms of Cm( z ) an d Cmm'( z ) may be found with 
the help of the explicit expression (p.21|) for Gq(z) and the integrals of 



Ref. [17|. We state only the explicit result for M(t), which is, for t — > oo, 
3^+i t r i . 

M (t) - ; fc+ i ™fc + i t (1 - 7)(* + 1 W + 3m ^+i 

log + ct 1 log ct V 



1 

log 2 ct 



+ 7— 2 —({l-±n 2 - 1 + ^ 2 )(k + l)(k + 2)m k 



3(1 - 7)(fc + 2)m k+1 + 9m fc+2 



(5.5) 



in which 7 = 0.577215... denotes Euler's constant. Both M(t)M'(t) and 
the product M(t) M'(t) the appear as t 2 times a power series in 1/ log ct of 
which the leading term is of order k + k' + 2, and in which the three leading 
orders have to be retained. Upon carrying out the subtraction one finds that 
the two leading orders cancel and the correlation AM(t)AM'(t) appears to 
be proportional to t 2 / log fc+fc +4 ct. Explicitly, as t — > 00, 



AM(t)AM'(t) ~ A 2 (k + l)(k' + l)rn k m' kf 3k+ k ^lf (t - 00; 

log + + ct 

in which the coefficient .4 is given by 



(5.6) 



A 2 = 1 + (F 1)3 - ±)vr 2 (5.7) 



Numerical evaluation gives Fi 3 = 0.27415..., whence^ 2 = 2.0608.... Eq. <K 
is the same as for the two-dimensional simple random walk |l4|, [?], ||, [i7j| , 
but with a different constant „4. 

6 Universality of fluctuations 

We consider in this section the normalized deviations from average 

M0 = ^L («) 



AM 2 (t) 
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These random functions of time satisfy by construction 



e M (t) = o, 



i 



(6.2) 



We consider now two arbitrary observables M and M'. When g < a < 1 
we have from Eq. ( |6,lD together with either Eq. ( |5.2[) or Eq. ( 3.13 ) that 



e M {t)e M '{t) 



i 



(6.3) 



It then follows from Eqs. fl6.2|) and (16^) that the difference 9 m — @M' is a 
random variable of zero average and zero variance. Such a random variable 
can only be itself equal to zero. We therefore deduce that, when § < a < 1, 
in the limit t — > oo all 0m{P) are equal to a single random variable, which 
we will call Q a (t), thus indicating explicitly its a dependence. 

When 1 < a < 2 we have for two observables M and M' whose orders, 
k and k' , are equal from Eq. (|6.l| ) and either Eq. (|5.6| ) or Eq. ( |4.10| ) again 
the result ( |6.3| ). Hence, when 1 < a < 2, in the limit t — > oo all 9m (t) with 
k = are equal to a single random variable - that we will call @ a o(t) -, and 
all 0M(t) with k = 1 are similarly equal to a single random variable - that 
we will call Q a i(t)- The variables @ a (t), Q a o(t), and @ a i(t) are universal 
in the sense that they are independent of the observables M (but depend at 
most on their order). 

In each of these cases the key ingredient necessary for arriving at Eq. 



is the factorization of M(t)M' '{€) into an M and an M' dependent part. This 
also explains why for a < | the same reasoning fails. 

The cross correlation between & a o(t) and ® a i(i) is easily found from the 
correlation between a 0M(t) and a 0M'(t) with k = and k' = 1, and use 
of the second one of Eqs. (gig). The answer is independent of the choice of 
M and M' , as it had to be, and reads 



©ao(*)@al(i) — # a lo(£>aOo£>all) 



(6.4) 



The coefficient ratio on the RHS of this equation depends only on the ex- 
ponent a and must necessarily be less than unity. In the limit a — ► 1 + it 
approaches unity and for a < 1 the distinction between surface and bulk 
observables is no longer reflected in the fluctuations. 

Upon combining all these conclusions we get explicitly 



AM(t) 



m 2 C(t\ogt)2 9|(i) 
m 2 B a t 2 -^ @ a (t) 
m k A(k + l)t(l log ct)- k - 2 0i (t) 

^ mi B a u t«~ l @ a l(t) 



(«=!) 
(| < a < 1) 
(a = l;k = 0,1) 

(Ka<2;fc = 0) 
(Ka<2;fc = 1) 



(6.5) 



Here all M dependence is contained in the coefficients m n . 
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7 Conclusions 



We have studied a large class of properties M(t) of the support of the one- 
dimensional t step Riemann walk. These include the number S(t) of distinct 
sites visited, and the number 1(f) of sequences of visited sites. The M{t) 
fall into two classes, the bulk or S-like properties, and the surface or I-like 
properties. The asymptotic laws found in the preceding sections for the 
averages, variances, and correlation of S(t) and I{t) have been summarized 
in Table I in the Introduction. 

It appears from that table that in the exponent regime < a < 1 

1/2 1/2 

the ratios AS 2 (t) / 5(f) and AI 2 (t) / /(f) tend to zero when t — ► oo, 
which indicates that the distributions of 5(f) and of 1(f) become infinitely 
narrowly peaked around their average. Hence in this exponent regime the 
ratio s(f) = 5(f)/I(f) represents the average number of sites per visited 
sequence. When a is strictly less than unity we have explicitly 

mi LSI G(0,1) -6(1,1) , 

lim alt) = = > 1 < a < 1 (7.1 

mii] G(0, 1) + G(l, 1) 



The first equality is based on Eq. (|3.5|) and in the second one we used the 



definition (3.4) of the m n [M] and the remarks at the end of Sec. 2.2. The 
finiteness of the result (|7.l| ) means that in the large t limit every new step 
of the walk creates a new visited sequence with a finite nonzero probability. 
This explains that in this exponent regime the asymptotic power laws do 
not distinguish between bulk and surface properties. For a — > 1~ expression 



(7.1) diverges, and when a = 1 the ratio s(t) increases logarithmically with 
t. 

In the exponent regime 1 < a < 2 the appropriately scaled distributions 
of S(t) and I(t) are of finite width even in the limit t — » oo. The support has 
an "interior", bulk and surface properties have different asymptotic power 
laws, and s(t) —* oo as t —* oo. In the terminology of critical phenomena, 
this regime is fluctuation dominated. In this regime the universality of 
fluctuations holds in a slightly weaker but at least as interesting a sense as 
for < a < 1. To describe the fluctuations, not a single but two universal 
stochastic variables are needed, one applying to the bulk and the other to 
the surface properties. These two variables become fully correlated in the 
limit a — > 1 + . 

Finally we remark that when 2/ a is equal to one of the integers 2, 3, 4, ... , 
the asymptotic laws of Table I coincide with the ones known to hold for the 
simple random walk on a lattice of dimension d = 2/ a. Similarly, the univer- 
sality properties for those a values have their analogs in the (i-dimensional 
simple random walk. Hence the "rule of the effective dimensionality", which 
states the correspondence a <3> 2/d, applies to all properties that we have 
studied. Of course it must break down when the comparison between the 
Riemann walk and the simple random walk is refined sufficiently. Also, we 
have not considered the borderline value a = 2, which is special ||, and for 
which this rule fails. 
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A Relations for the inverse sums and 

We collect here some elementary matrix algebra relations useful for dealing 
with the inverse sums Ga and qa occurring in the main text. The z depen- 
dence of these quantities plays no role. The presentation and notation are 
independent of the body of the paper. 

Let L be an invertible £ x I matrix. We define the "inverse sum" 1{L) 

by 

1-\L) = Y J ^ (A.l) 
In the remainder a, (3, and 7 will denote constants. 

Property 1. Let J be the £ x £ matrix with all Jij = 1, and let M be an 
invertible £ x £ matrix. Let L = aJ + jM. Then 

J(L) = a + tZ(M) (A.2) 

The proof of this relation is given in Ref. [|l7| . 

Property 2. Let M and N be invertible matrices of dimensions m x m 
and n x n, respectively, and let L be the block diagonal £ x £ matrix with 
blocks M and N. Then 

111 X 
+ T7at\ ( A -3) 



1(L) T(M) 1(N) 



This follows directly from the definition ( |A.1| ). The calculation of I(L) for 
an £ x £ matrix may be reduced to an inversion problem of dimension less 
than £ also in certain cases where L is not block diagonal, as shown below. 

Property 3. Let J mn be the m x n matrix with all elements equal to 1. 
Let L be £ x £ and of the form 

T ( 7 M (3J mn \ 



Then 



7 2 1(M)1(N) - p 
1[L) ~ 1 l{M)+ 1 l(N)-2(3 {A - b) 
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To prove this we rewrite L as L = (3 J + L, where J is as before and where 

~ / 7 M-pj mm \ 

L = [ 7 N - f3J nn ) (A - 6) 

From Property 1 we have that 2(L) = (3 +I(L), after which by applying 
Property 2 and once more Property 1, we obtain after some rearrange- 
ment Eq. ( |A.5f) . For = Eq. ( |A.5|) reduces to Property 2. 

In this work the need for PROPERTIES 1 and 3 arises when the limit 
7 — > has to be taken. For 7 = the matrices L that occur on the LHS 
of Eqs. (|A^D and Q are no longer invertible, but these properties allow 
nevertheless X(L) to be calculated in that limit. 



B Corrections to scaling 



In Eq. (|3.6| ) we have neglected the 0(V 2 ) terms that appear in Eq. ( 2.15| ) 



Since in the last step that led to Eq. ( |3.13| ) the leading order in 1 — z went 
down due to cancellations, we must now check that the 0(V 2 ) terms remain 
subdominant. In this Appendix we will write Eq. ( |2.13| ) in the simplified 
notation G^y^+s) = G + W where G and W are the first and second 
matrix, respectively, on the RHS of Eq. ( p. 13 ). Upon writing the inverse 



C*Au(r+B) as a perturbation series in W and applying Eq. ( |2.4D one finds 

00 

Gau(, + b) = £(-!)'£ [G-^WG-YW (B.1) 

e=o c ,c' 

The t = term of this series is the term shown explicitly on the RHS of 
Eq. ( |2.15| ), and has been the object of study in Sec. |3.2| . We will show here 
that the terms with I > 1 produce, in the scaling limit, only higher order 
corrections to the final result. To this end we first consider Cmm'(z) defined 
by Eq. fl2.6|) . Let Re(z) denote the contribution to Cmm'{ z ) from the £th 
term in Eq. ( |B.l| ). In order to estimate the order in 1 — z of Rg{z) as 
z — > 1 we first deduce from Eqs. ( 2.14| ) and ( [2.25|) that in the scaling limit 



the matrix elements of W behave as V at r+b — (1 — z) <* l (b — a) F'(£), and 
that summing on r amounts to applying (1 — z)~a j d£. This yields the 
asymptotic proportionality 

where Gq(z) represents the order in 1 — z of the matrix G. When I is odd, 
this integral vanishes by symmetry, which shows that the leading correction 



is of order V 2 , as anticipated. For £ — > we have, in virtue of Eq. ( 2.27 ) 



that F'(£) ~ £ 2+a . Hence the £ integral in Eq. (B.2) diverges in the origin 



for all £ > 1 when < a < 1. This signals that the main contribution 
comes from r values near the origin. The order in 1 — z of Ri(z) may then 
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be estimated by introducing in the integral the cutoff |£| ~ est x (1 — z)<*, 
which leads to Re(z) ~ (1 — z)°. When £ > 1 these additive corrections to 

2 o 

Cmm'(z) in Eq. (3.11) are negligible, therefore, with respect to the (1— z)~~ 
term which, in the relevant exponent regime | < a < 1, determines the final 
result. 
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